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Spectral singularities are predicted to occur in a non-Hermitian extension of the Friedrichs-Fano- 
Anderson model describing the decay of a discrete state \a) coupled to a continuum of modes. A 
physical realization of the model, based on electronic or photonic transport in a semi-infinite tight- 
binding lattice with an imaginary impurity site at the lattice boundary, is proposed. The occurrence 
of the spectral singularities is shown to correspond either to a diverging reflection probability (for 
an amplifying impurity) or to a vanishing reflection probability (for an absorbing impurity) from 
the lattice boundary. In the former case, the spectral singularity of the resolvent is also responsible 
for the non-decay of state \a) into the continuum, in spite of the absence of bound states. 
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INTRODUCTION. 



Over the last decade, a great attention has been de- 
voted to investigate the properties of non-Hermitian 
physical systems. In particular, it has been shown that 
the framework of quantum mechanics can be extended 
by relaxing the common constraint of Hermiticity for 
the underlying Hamiltonian H (see, for instance, [H-Il] 
and references therein), provided that H has a real en- 
ergy spectrum and is diagonalizable. In this case, af- 
ter a proper change of the inner product (metric) of the 
Hilbert space, the non-Hermitian Hamiltonian H may 
be used to define a unitary quantum system Exam- 
ples of such systems include non-Hermitian Hamiltonians 
with parity-time (VT) symmetry in which the real- 
ity of the energy spectrum (bound as well as radiation 
states) below a symmetry-breaking transition has been 
proved for several complex potentials [1]. Unfortunately, 
non-Hermitian Hamiltonians possessing a real- valued en- 
ergy spectrum may fail to be diagonalizable because of 
the occurrence of exceptional points in the point spec- 
trum 0, Q , or of spectral singularities in the continuous 
part of the energy spectrum 0-Il|- Exceptional points 
refer to the coalescence of two or more discrete eigenval- 
ues together with their eigenvectors; their physical rele- 
vance has been investigated in several works, and differ- 
ent physical realizations have been proposed and experi- 
mentally demonstrated JUj . On the other hand, spectral 
singularities refer to divergences of the resolvent operator 
G{z) = [z — H)~^ belonging to the continuous spectrum 
of H , i.e. which do not correspond to square-integrable 
eigenfunctions. As opposed to exceptional points, spec- 
tral singularities have received less attention from physi- 
cists and have been mainly viewed as a curious mathe- 
matical property of certain non-Hermitian operators !7r 
i]. The physical meaning and relevance of spectral sin- 
gularities have been highlighted solely quite recently no- 
tably by A. Mostafazadehlll] (see also H, i, [H, [ij) in 
the framework of wave scattering by complex potentials 
Mostafazadeh showed that spectral singularities of 
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a non-Hermitian Hamiltonian with a complex potential 



correspond to divergences of reflection and transmission 
coefficients of scattered states, i.e. to resonances with 
vanishing spectral width. He also investigated in details 
the appearance of spectral singularities in an electromag- 
netic realization of a non-Hermitian Hamiltonian with 
VT symmetry based on a waveguide filled by an atomic 
gas [Tllj [3 , following an earlier proposal by Ruschhaupt 
and coworkers [l^. 

It is the aim of this work to investigate the onset of 
spectral singularities in a non-Hermitian extension of 
the famous Friedrichs-Fano-Anderson (FFA) model [TtI - 
[Tgj . which generally describes the decay of a discrete 
state coupled to a continuum of modes. The FFA 
model is encountered in different areas of physics, rang- 
ing from atomic physics [20 22] to quantum electrody- 
namics [1^ HI] and condensed- matter physics [25l - [29| . 
The FFA model has been studied in the quantum the- 
ory of non-integrable systems [s^l and used to describe 
unstable quantum systems, quantum mechanical decay 
and quantum Zeno dynamics (sil . l32j . Simple models of 
single-particle electronic or photonic transport in tight- 
binding lattices can be also described by means of FFA 
Hamiltonians |25| - |29| . [33k38j . In this work it is shown 
in particular that spectral singularities can be observed 
in photonic or electronic transport in semi-infinite tight- 
binding lattices with an imaginary impurity site at the 
lattice boundary, leading to either a diverging or a van- 
ishing wave reflection from the lattice boundary. 
The paper is organized as follows. In Sec. II we intro- 
duce a non-Hermitian extension of the FFA model, and 
derive the conditions for the appearance of spectral sin- 
gularities. In particular, it is shown by direct calcula- 
tions that the appearance of divergences of the resolvent 
operator G{z) on the branch cut correspond to the non- 
diagonalizability of the FFA Hamiltonian. In Sec. HI we 
present an example of non-Hermitian FFA model show- 
ing spectral singularities, which provides a simple model 
of electron or photonic transport in a semi-infinite tight- 
binding lattice with a boundary impurity site. It is shown 
that spectral singularities in this model correspond to ei- 
ther a vanishing or a diverging reflection probability from 
the lattice boundary. Finally, in Sec. IV the main conclu- 
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sions are outlined. 



II. 



SPECTRAL SINGULARITIES IN A 
NON-HERMITIAN 
FRIEDRICHS-FANO- ANDERSON MODEL 

A. The model 

The standard FFA model describes the interaction of 
a discrete state \a), of energy Ea, with a continuous 
set of states |A:) with energy E{k) (see, for instance, 
[H [H-ii, m, [si]). Here \a) and \k) represent a com- 
plete set of Dirac states in the Hilbert space which sat- 
isfy the orthonormal conditions (a|a) = 1, (a|fc) = and 
(fc'jfc) — 6{k — k'). The Hamiltonian of the full system 
can be written a.s H = + V, where 



Ho = Ea\a){a\ + / dk E(k)\k){k\ 



(1) 



is the Hamiltonian of the non-interacting discrete and 
continuous states, and 



V= dk[v{k)\a){k\+v*{k)\k){a\] 



(2) 



is the (self-adjoint) interaction term, described by the 
spectral coupling function vlk). Typically, we assume 
that the energy E{k) of continuous states spans the in- 
terval El < E < E2 (eventually E2 = 00), and E{k) is 
a monotonic function (either increasing or decreasing) of 
A;, i.e. we assume that there are not energy degenera- 
cies of continuous states. The state vector of the system 
lip) evolves according to the Schrodinger equation (with 



' dt 



m 



(3) 



Note that, provided that the energies Ea and E{k) are 
real-valued, the Hamiltonian H is Hermitian. The spec- 
trum of H can be determined by either analyzing the 
singularities of the resolvent G{z) = {z — H)~^, or by 
projecting the eigenvalue equation i?|f/') = on the 

basis {|a), |fc)}. As is well known, the continuous spec- 
trum of if is -El < E < E2, i.e. the same as that of Hq, 
whereas the point spectrum can be either empty or com- 
posed by a number of discrete eigenvalues, either outside 
or embedded into the continuous spectrum (see, for in- 
stance, Hi,!!!,!!^). 

We now relax the Hermiticity condition of the FFA model 
by allowing the 'energy' Ea of the discrete state \a) to 
be complex valued. However, we will assume that the 
spectrum of H remains real- valued in spite of the non- 
Hermiticity of H. The condition for the spectrum of H 
to remain real-valued will be discussed below in Sec.H.B 
and corresponds to the absence of bound states, i.e. to an 
empty point spectrum. This means that a non- Hermitian 



FFA Hamiltonian has a real-valued energy spectrum if 
and only if its spectrum is purely continuous. We are 
interested here to determine, if any, the appearance of 
spectral singularities of H, which would prevent H to be 
diagonalizable. This problem can be addressed in two 
ways: (i) by the determination of the resolvent G{z), 
and (ii) by a direct diagonalization of H following the 
original procedure by Fano [18], extended to account for 
the non-Hermitian nature of H. In the former case, a 
spectral singularity at £ — Sq, embedded in the continu- 
ous spectrum (i?i,i?2), is revealed as a divergence of the 
Green function Gix,y;z) = {x\G{z)y) in a neighborhood 
of z — £q, divergence which does not corresponds to a 
bound state embedded in the continuum Q. In the lat- 
ter case, a spectral singularity ed £ — Sq occurs when 
{£q\£o) — 0, where |£^o) and |£^g) are the eigenf unctions 
of H and of its adjoint H\ respectively, corresponding 
to the eigenvalue £0 0. 



B. Spectral singularities: the resolvent approach 

For a given complex number z, the resolvent operator 
G{z) of the Hamiltonian H is defined as 



G{z) = {z-HY 



(4) 



i.e. G{z){z — H) = [z — H)G{z) =1, where I is the iden- 
tity operator. The knowledge of the resolvent of H, for 
any z, is equivalent to the knowledge of the set of eigen- 
functions and eigenvalues of H. In particular, the singu- 
larities of G{z) in the complex plane define the spectrum 
of H: an eigenvalue £ belonging to the point spectrum of 
H is a pole of G{z), whereas the branch cut of G{z) de- 
termines the continuous part of the spectrum of H. The 
Hamiltonian H is said to have a spectral singularity at 
£ = £q, where £q belongs to the continuous spectrum of 
H , if the function 



(5) 



is unbounded in the neighborhood of z = £q, and £0 
does not belong to the point spectrum of H, i.e. it does 
not corresponds to a bound state embedded in the con- 
tinuum In the previous equation, |x) and \ip) are 
two assigned functions of the Hilbert space; in particular, 
for \x) = \x) and \ip) = |y) one obtains the coordinate 
representation of the resolvent G, i.e. the Green func- 
tion Q{x,y;z) — {x\Giz)y). An interesting property of 
the FFA Hamiltonian is the possibility to calculate the 
resolvent in a closed form. The procedure to calculate 
G{z) is well known for the Hermitian case (see, for in- 
stance, (26j). and can be extended mutatis mutandis to 
the non-Hermitian FFA model considered in this work. 
As detailed in the Appendix A, the matrix elements of 
the resolvent G{z) on the complete basis {|a), read 
explicitly 
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g^^^ = {a\Giz)a) ■■ 

g,^k = {a\G{z)k) ■■ 
gk,a = {k\G{z)a) 

gk,k' = {k\G{z)k') ■■ 

where E(2;) is the self-energy, defined by 



1 



z-Ea- S(z) 



{z-E{k)){z-E^-i:{z)) 

v*{k) 

{z-E{k)){z-Ea-^{z)) 
v{k')v*{k) 



+ 



5{k - k') 

{z- E{k)){z- E{k')){z- Ea-E{z)) ' z-E{k') 
I 



dk 



z- E{ky 



(10) 



Note that, by introducing the density of states p{E) = 
{dE/dk)-^ and letting V{E) = p{E)\v{E)\^ {V{E)=0 
for E > E2 and E < Ex), the self-energy can be written 
in the equivalent form 



-f 

J El 



E 



(11) 



Note also that is not defined on the segment 

{Ei,E2) of the real axis, and one has 



E{z = £± iO+) = A{£) T mV{£) 



(12) 



where 



A(£:) = P 



E2 



dE 



El 



VjE) 
£-E' 



(6) 
(7) 
(8) 
(9) 



(13) 



P denotes the principal value, and £ is real-valued. 
Therefore, S(2) has a branch cut in the interval {Ei,E2)- 

For two assigned functions |x) = Xa|a) + / dkx{k)\k) 
and |(^) — ipa\a) + J dkip{k)\k) of the Hilbert space, the 
complex function gx,ip{z) = {x\G{z)(p) can be readily 
calculated as 



GxAz) = XlVagaAz) + X*aj dkip{k)ga,k{z) + V>a j dkx* {k)gk,a{z) + J dkdk' X* {k')^{k)gk' ,k{z) (14) 



I 

Substitution of Eqs.(6-9) into Eq.(14) finally yields where we have set 

gxAz) = Ga,a{z)^l{z) + ^2{z) (15) 
I 



^ ^1 ,^p{E)v(EME)^ ^1 ^^p{E)v*iE)r{E)^ ^^^^ 

Mz) = [ dE^-i^)2<*i^. (17) 
J z — E 



We are now ready to determine the singularities of the singularities. To this aim, let us notice that ^i{z) and 
resolvent, i.e. the spectrum of H and possible spectral $2 (-2) are bounded functions of z and have a branch cut 
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on the segment {Ei,E2) of the real axis (as for the self- 
energy E). According to Eq.(15), we may therefore Umit 
to consider the singularities of Ga,a{z)- From Eqs.(6), 
(11) and (12) one can conclude that: 

(i) The continuous spectrum of H is the same as that of 
i/o, i-e. the interval {Ei,E2) of the real axis, where G{z) 
has a branch cut. 

(ii) The point spectrum of H, corresponding to bound 
states outside the continuum, are the complex roots z of 
the equation 



Z^Ea^ S(z) 



(18) 



in correspondence of which the resolvent G{z) has a pole, 
(iii) A spectral singularity £q of the continuous 
spectrum(ii'i, £^2) is any solution of the coupled equa- 
tions 



lm{Ea) = ±ttV{£„) 
£o-Re{Ea) = A(fo) 



(19) 
(20) 



provided that lm{Ea) 7^ [4l|. 

Note that real- valued energies can not belong to the point 
spectrum of H because Eq.(18) does not have real-valued 
roots whenever H is non-Hermitian [l^ . This means that 
the point spectrum of H, if not empty, is strictly complex- 
valued. On the other hand, the continuous part of the 
spectrum is real- valued according to the property (i). 
Therefore, we may conclude that the energy spectrum 
of the non-Hermitian FFA Hamiltonian is real-valued if 
and only if its spectrum is purely continuous. In this 
case, spectral singularities in the continuous spectrum 
occur whenever Eqs.(19) and (20) can be simultaneously 
satisfied. It should be noted that the behavior of G{z) 
in the neighborhood of the spectral singularity z = Sq 
on the continuous spectrum is different for an 'absorb- 
ing' [Im(£;a) < 0] and for an 'amplifying' [Im(i?a) > 0] 
complex energy Ea of state \a). Since the spectral cou- 
pling V{£o) is always positive, for an absorbing complex 
energy Eq.(19) can be satisfied by taking the lower (neg- 
ative) sign on the right hand side; correspondingly, from 
Eqs.(6) and (12) it follows that Ga.a{z) is unbounded 
when z ^ £q from the bottom of the real energy axis, 
i.e. for Im(2:) < 0, but Ga,aiz) remains bounded when 
z ^ £0 with Im(z) > Conversely, for an 'amplify- 
ing' complex energy lm{Ea) > 0, Eq.(19) can be satisfied 
by taking the upper (positive) sign on the right hand side; 
correspondingly, Ga,a{z) is unbounded when z — ^ £0 with 
Im(z) > 0, but it remains bounded when z — >■ with 
Im(z) < [4^. Such a different behavior of spectral 
singularities for an absorbing or an amplifying complex 
energy Ea has some relevant physical implications, that 
will be discussed in Sec. II. D and Sec. III. C. In Sec. II. D it 
will be shown that the appearance of a spectral singular- 
ity in the amplifying case is responsible for the non-decay 
of state I a) into the continuum, in spite of the absence 
of bound states; in Sec. III. C the interplay between spec- 
tral singularities and wave scattering will be investigated 
for a semi-infinite tight-binding lattice realization of the 
FFA Hamiltonian. 



C. Spectral singularities: the Fano diagonalization 
procedure 

Let us assume that the spectrum of the non-Hermitian 
Hamiltonian H = Hq -|- is real- valued. As shown in 
the previous subsection, this implies that the spectrum 
of H is purely continuous and spans the interval {Ei, E2). 
Let us indicate by \£) the (improper) eigenfunction of H 
corresponding to the eigenvalue £, and by \£'^) the (im- 
proper) eigenfunction of the adjoint corresponding 
to the same eigenvalue £. For H to be diagonalizable, 
the set of functions {|f ), \£^)} must form a complete 
biorthonormal basis of Hilbert space , that is 

{£^\£') = 5{£ - £') , [ \£\£){£^\^I. (21) 

JEi 

A spectral singularity at £ ~ £0 sets in when 



{£o\£l)=o, 



(22) 



which prevents H to be diagonalizable. To determine the 
onset of spectral singularities of the FFA Hamiltonian H, 
we can thus proceed by calculating the eigenfunctions \£) 
of H and jf^) of its adjoint H'' following the diagonal- 
ization procedure used by Fano [l^ in the problem of 
atomic autoionizing resonances (see also [13, El), prop- 
erly modified to account for the non-Hermitian nature of 
H. To this aim, we expand the eigenstates \£) and \£'') 
on the orthonormal and complete basis {|a), \k)} as 



\£) = a{£)\a)+ J dk(3{£,k)\k) (23) 
= a^£)\a)+ j dkp\£,k)\k) (24) 



with expansion coefficients a{£), a^{£)^ fi{£,k) and 
/3''{£,k) to be determined. Note that 



{£'\£^) = a*i£')a^£) + J dkp*{£', k)P\£, k). (25) 

Using Eqs.(l), (2) and (23), the eigenvalue equation 
H\£) = £\£) yields the following coupled equations for 
the expansion coefficients a{£) and /3{£,k): 

{Ea - £)a{£) + j dkv{k)p{£, k) = (26) 

[E(k)-£](3{£,k)+v*{k)a{£)^0. (27) 

Equation (27) can be solved for (3{£,k) and substituted 
into Eq.(26). However, there is a singularity at E{k) — £, 
so that 1/ [E{k) — £] must be written as its principal and 
delta function parts, yielding (see [isl [21}) 



P{£,k) = -v*ik)a{£) P 



1 



E{k) - £ 



z{£)5{E{k) - £\ 



(28) 

The coefficient z{£) in front of the delta function on the 
right hand side of Eq.(28) is determined by substituting 
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Eq.(28) into Eq.(26) and requiring that a{£) does not 
vanish. This yields 



zi£) 



Ea^£ + A(g) 
V{£) 



(29) 



where A{£) is defined by Eq.(13), ¥{£) = p{£)\v{£)\^ 
and p{£) — {dE{k)/dk)^^ is the density of states. 
Similarly, the eigenvalue equation H^£^) = £\£'') yields 



{E*-£)a^£)+ J dkv{k)l3''{£,k)^0 (30) 
[E{k) ~ £] I3\£, k) + v*{k)a'^{£) = (31) 

and the following expression of the coefficient /?t(£,fc) 
can be derived following the same procedure as above 

/3t(£-, k) = -v*{k)aH£) (P ^(^|_g + z*{£)5{E{k) ~ £) 

(32) 

Substitution of Eqs.(28) and (32) into Eq.(25) yields 

{£' \£^ )= a* {£')a\£)F {£,£') (33) 
where we have set 

F{£,£') = dEV{E)(p^^ + z*{£)S{E-£)^ 

X (p^^ + z*i£')S{E-£')^+l. (34) 

The calculation of the integral on the right hand side of 
Eq.(34) is complicated by the product of two principal 
parts, which must be properly resolved into its partial- 
fraction and singular terms (see [isj or the Appendix of 
Ref.[2l[). Taking into account that 



P^P ' 



= P- 



1 



1 



-P- 



1 



E-£ E-£' £-£'\ £'-E £-E^ 

+ TTHi£' - E)6{£ - E) (35) 

and using Eq.(29), from Eq.(34) one obtains 

F{£,£') = V{£) [tt^ + z*^{£)] 5{£ - £') (36) 

so that [see Eq.(33)] 

{£'\£^) =a*{£)a\£)V(£) [tt'^ + z*^{£)\ 5{£-£'). (37) 

In order to represent a complete biorthonor- 

mal set of functions [Eq.(21)], the amplitudes a[£) and 
a^(£) should be thus normalized such that 



a*{£)a\£)V{£) [n'^ + z*''{£)\ = 1. 



(38) 



For the Hermitian FFA model, the energy Ea is real- 
valued, z{£) turns out to be real-valued and a^{£) — 
a{£) [because of \£'^) — [£)], so that Eq.(38) is solved to 
yield 



V{£) 



n^V^{£) + [Ea-£ + A{£W 



(39) 



according to Fano [l^. The physical meaning of Eq.(39) 
is well known [l8l |: owing to the coupling with the con- 
tinuum, the discrete state \a) is 'diluited' throughout 
a set of continuous states (i.e., it becomes a resonance 
for H) with a typical resonance curve |q;(£)|^ peaked at 
£^Ea- A{Ea) of width ~ nV{Ea) M- 
For the non-Hermitian FFA model, the energy Ea is 
complex- valued and from Eq.(37) it follows that a spec- 
tral singularity can appear at the energy £ = £o oi the 
continuous spectrum such that tt^ -|- z*^(£'o) = 0, i.e. 
when 



z*{£o) 



(40) 



Using Eq.(29), Eq.(40) yields the following conditions 
for the appearance of a spectral singularity in the non- 
Hermitian FFA model 

Im(Sa) - ±^V{£a) , Re(£;,) -£a + A{£a) - (41) 

which are precisely Eqs.(19) and (20) derived in Sec. II. B 
following the resolvent approach to spectral singularities. 



D. Spectral singularities and decay dynamics 

Hermitian FFA models are often used to describe the 
evolution of unstable quantum systems and the related 
problem of quantum mechanical decay and deviations 
from an exponential decay law (see, for instance, [sil l32j | 
and references therein). If the system is initially prepared 
in state |a), i.e. if \'ip{t — 0)) — |a), the survival proba- 
bility P{t) to find the system at state \a) at a successive 
time t is given by P{t) = |ca(i)|^, where 



Ca{t) = {a\ij{t)) = (a|cxp(-zi/t)|a) 



(42) 



Here we consider the case of a non-Hermitian FFA Hamil- 
tonian H with a real- valued (and thus purely continuous) 
energy spectrum, and briefly discuss the physical implica- 
tions of spectral singularities of H on the decay dynamics 
of survival probability P{t). The temporal evolution op- 
erator exjp{—iHt) of the system can be written in terms 
of the resolvent G{z) as 



exp{—iHt) — — dzG{z) exp{—izt) 
27r Jb 



(43) 



where the Bromwich path B is any horizontal line 
Im(z) = const > in the upper half of the complex z 
plane [see Fig.l(a)]. Substitution of Eq.(43) into Eq.(42) 
yields 



Ca(0 = 77- / rfz^/a,a(^)cxp(-izi) 

27r Jb 



(44) 



where the matrix element Ga,aiz) of the resolvent is given 
by Eq.(6). As discussed in Sec. II. B, Ga,a{z) is analytic in 
the full complex plane, expect for the branch cut on the 
real axis, from z = Ei to z = E2 (see Fig.l), where it 
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(a) 



lm(z) 



(b) 



lm(z) 



^2 Re^) 



Re(z) 



FIG. 1; (a) Integration contour (Bromwich path B) in the 
lm(z) > complex plain entering in Eqs.(43) and (44). The 
bold horizontal segment on the real axis is the continuous 
spectrum of H and corresponds to a branch cut of Qa,a{z). 
(b) Integration contour (Hankel paths hi and /12) after de- 
formation of the Bromwich path. The integration along the 
solid (dashed) curves of the Hankel paths is made on the first 
(second) Riemannian sheet of Qa,a{z). 



can also become unbounded at energies £q corresponding 
to spectral singularities. The integral on the right hand 
side of Eq.(44) can be evaluated by the residue method 
after suitably closing the Bromwich path B with a con- 
tour in the Im(z) < half-plane (see, e.g., [31]). Since 
the closure crosses the branch cut on the real axis, the 
contour must necessarily pass into the second Rieman- 
nian sheet in the sector of the Im(z) < half-plane with 
El < Re(z) < E2, whereas it remains in the first Rie- 
mannian sheet in the other two sectors Re(2:) < Ei and 
Re(2;) > E2 of the Im(z) < half-plane. To properly 
close the contour, it is thus necessary to go back and 
turn around the two branch points of the cut at z = Ei 
and z = E2, following the Hankel paths hi and /12 as 
shown in Fig. 1(b). Indicating by Qa^P{z) the analytic 
continuation of Ga,a{z) from the upper [Im(2) > 0] to 
the lower [Im(z) < 0] half-plane across the branch cut, 

and by zi, Z2, ■■■ the (possible) poles of Gi/,a\z) that hes 
in the sector Im(z) < ij, Ei < Re{z) < E2 of the complex 
plane (with rj > Q arbitrarily small), assuming that the 
poles are of first order one can write 



bounded function. Therefore P{t) ^ as < — ^ cx) t45|. 
Typ ically, Q^.a\z) has one pole below the imaginary axis 
[3l|, so that the decay of Ca{t) follows an exponential law 
(with decay rate determined by the imaginary part of 
the pole), corrected by the cut contribution C{t). Let 
us consider now the non-Hermitian FFA model with a 
spectral singularity a,t £ = £q in the continuous spec- 
trum. In this case, we have to distinguish two cases. For 
an absorbing complex energy Ea^ as shown in Sec. II. B 
Qa,a{z) is bounded for z — E + iQ'^\ therefore, its analytic 

continuation Qa,a' [z] in the lower-half complex plane is 
bounded for z = — iO+. This means that possible 
poles of ^a^i"* {z) have a strictly negative imaginary part, 
as for an Hermitian FFA model. Therefore the survival 
probability P{t) decays toward zero as i — > 00 similarly 
to what happens in the Hermitian limit. Conversely, for 
an amplifying complex energy Ea^ Qa,a{z) diverges for 

z = £q + iO+; therefore Q'i^^i\z) has a pole on the real 
axis, zi — £o- In this case, according to Eq.(45) the 
survival probability P(t) does not decay, and may be- 
come larger than one depending on the modulus of the 
residue of the dominant pole. This phenomenon is anal- 
ogous to the fractional decay found in Hermitian FFA 
models and related to the existence of bound states (i.e. 
a non-empty point spectrum of H). However, in the non- 
Hermitian FFA model with Ini(£'a) > the non-decaying 
behavior of P{t) results from the appearance of a spec- 
tral singularity in the continuous spectrum, not from the 
existence of bound states. It should be noted that such 
a non-decaying behavior was previously predicted in the 
study of the lasing threshold of an optical microcavity 
resonantly coupled to a coupled resonator optical waveg- 
uide under special coupling conditions [s^ . This unusual 
behavior of laser phase transition was explained as a non- 
markovian effect arising from the structured continuum 
of the decay channel, however it was not related to spec- 
tral singularities of the underlying Hamiltonian. 



Ca(i) ^Y.^f' exp(-izfct) +C{t) 



(45) 



where Rk is the residue of Ga^.P{z) at z = Zk, and C{t) is 
the contribution from the contour integration along the 
Hankel paths hi and /i2 

C{t) = — / dzlGaAz) - G^'Piz)] expi~izt) + 



2^ 



-E2 + ffl 



-E2 — too 



dz{g^^^X^)-QaA^)\^M-^zi\^^) 



The cut contribution C(t) vanishes as i — >■ 00, however 
it is responsible for the appearance of nonexponential 
features in the decay dynamics [3l| . For the Hermitian 

FFA model, the poles z^ of Q^^a^P {z) lie below the real 
axis, i.e. Im(zfc) < 0, because Qa.a[,z = E + iO~^) is a 



III. SPECTRAL SINGULARITIES IN A 
SEMI-INFINTE TIGHT-BINDING LATTICE 
WITH A BOUNDARY IMPURITY SITE 



In this section we present a simple and analytically- 
solvable example of a non-Hermitian FFA model showing 
spectral singularities, which describes rather generally 
single-particle electron or photon transport on a semi- 
infinite one-dimensional tight-binding lattice with an im- 
purity site. The model is first presented in the framework 
of the general theory developed in Sec. II, and its tight- 
binding lattice realization is subsequently described. The 
physical implications of spectral singularities on wave 
scattering from the lattice boundary and on the decay 
dynamics of the impurity site are finally highlighted. 
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A. The Hamiltonian 

Let us consider the non-Hermitian FFA model defined 
by the following relations for the energy dispersion E{k) 
and spectral coupling v{k) 



E{k) — — 2KoCOsfc , v{k) 



K„ sin k 



(47) 



where kq, Ka a-re two real- valued positive constants and 
< fc < TT. The Hermitian limit of this model, at- 
tained by assuming Im(i?a) = 0, is a special case of 
the FFA model previously investigated in Ref. [sB] , which 
is exactly solvable (see also [1^). Note that the con- 
tinuous spectrum of H spans the band (i?i,i?2), with 
E2 = —El = 2ko. The density of states for this model is 
given by 



PiE) 



dE 
'dk 




-2ko < E < 2ko 
\E\ > 2ko 



(48) 

which shows van-Hove singularities at the band edges, 
whereas the positive spectral function V{E), defined by 
V{E) = p{E)\v{E)\^, reads 



V{E) = 



'2ko < E <2kq ^4g^ 

\E\ > 2ko 




-1 1 
lnn(£-„)/Ko 



1 1 
lm(£j/Ko 



FIG. 2: Domains of non-existence of bound states for the 
Hamiltonian H in the {lm{Ea)/ t^o, i^a/ 1^0) plane (shaded re- 
gions) for increasing values of the ratio |Re(i5a)|/Ko- For a 
non-Hermitian Hamiltonian, i.e. lm{Ea) 7^ 0, in the shaded 
regions the energy spectrum of H is real-valued and purely 
continuous. Spectral singularities occur at the boundary of 
the shaded regions. 



which is non-singular. Substitution of Eq.(49) into 
Eg. (11) yields the following expression for the self-energy 

s(z) M 



S(z) = 



+ iz 



(50) 



and thus [see Eq.(12)] 
A(£) = Re(S(z = £±iO+)) 
4(^ + v/^ 



2k; 

h 

2 



(51) 

£ < -2kq 

— 2ko i£ £ i£ 2ko 

£ > 2ko 



The condition for the non-Hermitian Hamiltonian to pos- 
sess a real- valued spectrum (i.e. to avoid complex- valued 
energies arising from bound states outside the contin- 
uum) is derived in Appendix B. Precisely, let ^1^2 be the 
two roots of the second-order algebraic equation 



Kq 



(52) 



Then the Hamiltonian H has a real- valued energy spec- 
trum if and only if \£,i,2\ < 1- Figure 2 shows the domain 
in the plane (lm.{Ea)/ Kq, Ka/ kq) where H has a purely 
continuous energy spectrum for a few increasing values 
of the ratio \Ke{Ea)/ kq\- The domain lies in the sector 



/"^a/i^o ^ V2 and shrinks toward lm{Ea) / kq = Kq/kq = 
as \Re{Ea) / Ko\ 2~ . For \Re{Ea) / ko\ < 2, bound states 
do exist for any value of Ka/i^o and Im(ii'a)/Ko. The 
wider domain is attained for Re(£'a) = 0. In particular, 
for Re(£'a) = and Ka/no = ^/2, from Eq.(52) it follows 
that H has a real- valued energy spectrum provided that 



2ko < lia{Ea) < 2k.q. 



(53) 



Let us now consider the occurrence of spectral singulari- 
ties. According to Eqs.(19) and (20) and using Eqs.(49) 
and (51), a spectral singularity at energy £ = £q, in- 
side the interval {—2ko,2kq), is found provided that the 
following two equations are simultaneously satisfied 



Im(i?,) = ±^1 

Kq 



£o_ 

2kq 



(54) 
(55) 



For arbitrarily given values of Ea, Ka and kq, the 
above conditions are generally not satisfied [nowhere for 
£0 in the range (— 2ko,2ko)], he- the non-Hermitian 
FFA Hamiltonian is generally diagonalizable. Spec- 
tral singularities appear solely when a constraint among 
Iie{Ea)/ Kq, lm(Ea)/Ko and Kq/kq is satisfied. Let us 
first assume Ka/KQ strictly smaller that ^/2. In this 
case, a single spectral singularity, at the energy £q = 



8 



\a) 



|1> |2> |3> \n) 



In the Wannier representation, one can readily show that 



1^0 '^0 

FIG. 3: (color online) Schematic of a semi-infinite one- 
dimensional tight-binding lattice attached to a boundary im- 
purity site I a). 



Re(_Ba)/(l — '«^/2kq) [see Eq.(55)], is found provided that 



W{Ea) 



(2ko - Kq/ko)^ 



(56) 



It can be readily shown that the condition (56) defines the 
boundary of the domains shown in Fig. 2, i.e. a spectral 
singularity appears when the boundary of existence of 
bound states is approached. The case Kq/ko = is 
somehow singular. From Eqs.(54) and (55), for Kq/'^o — 
it follows that there are two spectral singularities at 
energies 



(57) 



provided that Re(i?a) = 0. The physical meaning of such 
spectral singularities will be discussed in Sec. III. C. 



B. Lattice realization 

FFA models are often encountered in connection to 
single-particle electronic or photonic transport in one- 
dimensional tight- binding lattices or networks (see, e.g., 
[25l-l2^. 's^-'ss'] and references therein) , and in most cases 
the underlying Hamiltonian is Hermitian. In particu- 
lar, the Hermitian limit of the FFA Hamiltonian H con- 
sidered in the previous subsection has been previously 
studied in [2^, [s^ and shown to be equivalent to a tight- 
binding Hamiltonian of a semi-infinity lattice with an 
impurity site. The equivalence can be proven after rep- 
resenting the Bloch states | k) of the tight-binding energy 
band in terms of localized Wannier states \n) on a lattice. 
Let us introduce the Wannier states 



\n) as 



\n) = \l — [ dksin{nk)\k) 



(58) 



(59) 





for n= 1,2, 3, .... Taking into account that 

dk sin(nfc) sin(mfc) = — (S„,m 
2 

(n, m > 1), one can readily show that the Wannier states 
form an orthonormal system, i.e. (n|m) — 5n,m- Addi- 
tionally, from Eq.(58) it follows that the Bloch states |A:) 
can be decomposed as a superposition of Wannier states 
|n) according to 



(60) 



\k) = \j— sin(nfc)|r 



dA:S(A:)|fc)(fc| =-Ato^ + 1| + \n+l){n\) 

n=l 

(61) 

and 

j dk[v{k)\a){k\+v*{k)\k){a\] = ^Ka{\a){l\ + \l){a\) 

(62) 

so that the Hamiltonian H = Hq + V can be written in 
the equivalent form 

oo 

H = -Kn5](l")(" + l| + l" + l)("l) + ^a|a)(a| + 



71=1 



-«:,(|a)(l| + |l)(a|). 



(63) 



In its present form, Eq.(63) describes single-particle elec- 
tron or photon transport on a one-dimensional semi- 
infinite tight-binding lattice [l^, [3^ |33| , with a hopping 
rate kq between adjacent sites of the lattice and with the 
boundary attached to an impurity site \a) with 'complex' 
potential energy Ea and with hopping rate Ka (see Fig. 3). 
From a physical viewpoint, the complex potential at the 
boundary impurity site may account for e.g. loss of the 
quantum particle flux into other decay channels (quan- 
tum absorbing potentials f??!), or optical gain or loss of 
light waves in photonic structures :34, 35j. For instance, 
light transport in a semi-infinite waveguide array as in 
Ref. 37], but with a lossy (or active) boundary waveg- 
uide, provides a simple and experimentally accessible re- 
alization of the tight-binding model (63). It should be 
noted that transport and scattering phenomena in tight- 
binding lattices with complex potentials have been the- 
oretically investigated in recent works (see, for instance, 
t47.,i48]), however spectral singularities were not found in 
these previous models. 



C. Spectral singularities and lattice wave reflection 

The physical meaning of spectral singularities in the 
non-Hermitian FFA model can be captured by analyzing 
the wave reflection properties of the lattice of Fig. 3. As 
it will be shown below, a spectral singularity corresponds 
to the appearance of a diverging peak in the reflectance 
spectrum when the boundary site is an 'amplifying' im- 
purity, i.e. when lm{Ea) > 0, and to the vanishing of 
wave reflection when the boundary site is an 'absorbing' 
impurity, i.e. when lm{Ea) < 0. In the former case, we 
retrieve for a 'discrete' scattering problem the physical 
explanation of spectral singularities as resonances with 
vanishing spectral width, established by Mostafazadeh 
Ref. 9] in the framework of 'continuous' wave scatter- 
ing by complex potentials. Conversely, the latter case, 
i.e. that of an 'absorbing' impurity site, gives a different 
manifestation of a spectral singularity: a wave incident 
on the lattice boundary is totally absorbed. 
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(b) 



10 



10 



10 
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Y 




3 



0.571 n 0.5 n iz 

wave number k 

FIG. 4: (color online) Probability of particle reflection R{k) 
for (a) an amplifying impurity site [Im(i5a) > 0], and (b) 
an absorbing impurity site [Im(i5a) < 0] for Ha/f^o = 1 and 
Re{Ea) = 0. Curves 1, 2 and 3 correspond to \Im{Ea)\/Ka = 
1/3, 2/3 and 1, respectively. According to Eqs. (54) and (55), a 
spectral singularity occurs at fe = 7r/2 (i.e., £o = — 2ko cos k = 
0) for \lTa{Ea)\/Ka = 1 (curve 3) 




0.571 71 0.57t 71 

wave number k 

FIG. 5: (color online) Same as Fig. 4, but for for Ka/Ko = 
y/2 and Ke{Ea) ~ 0. Curves 1, 2 and 3 correspond to 
\lm{Ea)\/ lia = 2/3, 4/3 and 2, respectively. According to 
Eq.(57), there are two spectral singularities, corresponding 
to the energies Sq = ±[4ko ~ Ini^(^'a)]^''^, which coalesce at 
|Im(£a)|/Ka = 2 (curve 3). 



To analyze the reflection properties of the lattice shown 
in Fig. 3, let us expand the state vector \ip{t)) as 



= ca{t)\a) +J2cn{t)\n) 



(64) 



n=l 



where |c„(i)p is the occupation probability of site \n) 
and |ca(t)p the occupation probability of the boundary 
impurity site \a). From Eqs. (3) and (63), it follows that 
the occupation amplitudes c„ and Ca satisfy the following 
coupled equations 



.dcn 
t—— 
dt 

.dci 
I—— 

dt 

.dCa 

I—— 

dt 



-Ko(Cn+l + Cn-l) 71 > 2 (65) 
-K0C2 - KaCa (66) 
-KaCi+EaCa- (67) 



Plane wave solutions to Eqs. (65-67) with wave number 
(momentum) fc (0 < /c < tt) , corresponding to eigenstates 
of H with energy E{k) = — 2KoCosfc, are of the form 
Cn{t) = Cn(k)eyiY>{-iE{k)t), Ca(t) = Cq(A:) exp(-ii5(fc)<), 
where 



Cn{k) 
Ca{k) 



exjp[—ik{n — 1)] 



— [exp(«fc) + r exp{~ik)] 



rik) exp[ik{n - 1)] (68) 
(69) 



and 



r{k) 



~ '*o(2 cos A: + Ea/nQj exp{ik) 
— Kq{2 cos k + Ea/Ko) exp(— ifc) 



(70) 



is the spectral reflection coefficient (see, for instance, 
[37|). Note that c„(t) is given by the superposition of 
two traveling waves, a regressive wave exp[—ikn — iE{k)t] 
that propagates along decreasing values of n, and a pro- 
gressive wave exp[ifcn — iE{k)t] that propagates in the 
opposite direction, i.e. which is reflected from the lat- 
tice boundary \4^. The probability of reflection from 
the lattice boundary is given by 



2kq) cos k 



: sinfc -I- Kolm{Ea)]'^ 



[{Kl~2Kl)cOSk- 



KoRe{EaW 
I 



: sinfc — Kolin{Ea)Y 



(71) 



In the Hermitian limit {lm{Ea) = 0), one has R{k) = 1, 
i.e. the incident wave is completely reflected from the lat- 
tice boundary. This is merely a consequence of conserva- 
tion of the particle probability in the scattering process. 
Conversely, for a complex-valued energy Ea of the impu- 
rity site, one has R{k) < 1 for lm{Ea) < 0, i.e. for an ab- 



sorbing potential, and R{k) > 1 for lm{Ea) > 0, i.e. for 
an amplifying potential. In particular, for Im(i?a) > 
the reflection probability R{k) goes to infinity at wave 
numbers fc = fco (with < ko < tt) that satisfy simulta- 
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neously the two conditions 

K^sinfco = KQlmiEa) (72) 
{nl~2nl)cosko = KoReiEa). (73) 

Physically, the condition i? — > cx) implies the existence 
of an outgoing wave in the lattice that is sustained by 
the amplifying complex potential at the impurity site. 
Such a divergence of R{k) is the signature of a spectral 
singularity of H. In fact, taking into account that the 
energy of incident/reflected waves is £q — ~2kq cos ko, it 
can be easily shown that Eq.(73) is equivalent to Eq.(55), 
whereas Eq.(72) is equivalent to Eq.(54) with the upper 
(positive) sign on the right hand side. Therefore, for 
an amplifying impurity site [lm{Ea) > 0], the condition 
i? — >■ oo is equivalent to the appearance of a spectral 
singularity. This equivalence extends, to our scattering 
problem on a truncated lattice, the general result shown 
by Mostafazadeh [9j , suggesting to interpret spectral sin- 
gularities of a non-Hermitian Hamiltonian as resonances 
with vanishing width. However, for an absorbing poten- 
tial energy at the impurity site, i.e. for Im(i?a) < 0, our 
lattice model indicates that the appearance of spectral 
singularities has a different physical interpretation. In 
fact, for Im(£'o) < the reflection probability R{k) is 
bounded from above and smaller than one, which pre- 
vents R{k) to diverge. However, in this case R{k) can 
vanish at wave numbers k = kg (with < ko < tt) that 
satisfy simultaneously the two conditions [see Eq.(71)] 

(74) 
(75) 

Note that Eq.(75) is equivalent to Eq.(55), whereas 
Eq.(74) is equivalent to Eq.(54) with the lower (nega- 
tive) sign on the right hand side. Therefore, for an ab- 
sorbing impurity site [Im(£'Q) < 0], the appearance of a 
spectral singularity is equivalent to the vanishing of the 
reflection probability R: an ingoing plane wave with wave 
number fco is fully absorbed by the impurity site at the 
lattice edge. An an example. Fig. 4 shows the behavior 
of R{k) for an amplifying [Fig. 4(a)] and for an absorbing 
[Fig. 4(b)] impurity site for Re{Ea) = and Ha/i^o = 1- 
The different curves in the figures refer to different values 
of Im(£'a)/Ko. Note that, at the value oilm{Ea)/ cor- 
responding to the appearance of the spectral singularity, 
a divergence and a zero in the R{k) curve are observed 
in Figs. 4(a) and 4(b), respectively. Figure 5 shows the 
behavior of R{k) as in Fig. 4, but for Re(i?a) = and 
i^a/i-^o = \/2. In this case there are two spectral singu- 
larities at energies given by Eq.(57), which explain the 
two peaks [Fig. 5(a)] or dips [Fig. 5(b)] in the reflectance 
curve R{k). By increasing \lva{Ea)\/ kq, the two peaks 
(or dips) get closer each other, until they coalesce at 
\lui{Ea)\/Ka - 2 (curve 3 

Figures 6 and 7 show two examples of wave packet 
reflection from the lattice boundary for an amplifying 
(Fig. 6) and an absorbing (Fig. 7) impurity site. The fig- 
ures show a snapshot of |c„(i)| as obtained by numer- 
ical analysis of Eqs. (65-67) assuming at i = a broad 



K^sinfco = -no^Yn{Ea) 
(k^ - 2Ko)cosfco = KoRe(£'a). 
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FIG. 6: (color online) Wave packet reflection in a semi-inflnite 
lattice with an amplifying impurity site for = Ka = 1, 
Re(£'a) = and \m{Ea) = 1. The initial wave packet is 
Gaussian- shaped, with peak amplitude |c„(0)| = 1 at n = 
no = 30, width An = 12 and momentum k — Stt/IO in (a), 
and k — TT /2 in (b) (corresponding to the spectral singularity 
energy £o — 0). In (b)the incident wave packet, indicated by 
an arrow, is not visible owing to the large amplification of the 
reflected wave. 
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FIG. 7: (color online) Same as Fig.6, but for an absorbing 
impurity site Im(_Ea) = —1 (other parameters as in Fig.6). 



Gaussian distribution of site occupation amplitudes, i.e. 
c„(0) = exp[— (n — no)^/An^ — ikn], where An is the 
wave packet width, k the mean wave packet momen- 
tum, and no ^ An the wave packet center of mass. 
The large (diverging) amplification of the reflected wave 
packet in Fig. 6(b), and the almost absence of wave packet 
reflection in Fig. 7(b), are clearly visible when the energy 
£ — — 2Kocosfc of the incoming wave packet attains the 
spectral singularity point £o = 0- 

The interplay between spectral singularities and de- 
cay dynamics, discussed in Sec. II. D, is exemplifled in 
Figs. 8 and 9. The system of Eqs. (65-67) has been nu- 
merically integrated with the initial condition Ca(0) = 1 
and c„(0) = 1. Note that, in the photonic realization 
of the semi-inflnite lattice model of Ref. [13] , such an ini- 
tial condition simply corresponds to initial excitation of 
the boundary waveguide. The behavior of the site oc- 
cupation probability Pit) — \cait)\'^ for an amplifying 
and for an absorbing impurity site is depicted in Figs. 8 
and 9, respectively. Note that, according to the gen- 
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t 

FIG. 8: (color online) Evolution of occupation probability 
P{t) = \ca{t)\'^ in a semi-infinite lattice with an amplifying 
impurity site for Re(_Ea) = 0, \m[Ea) ~ 1 and for (a) kq — 

Ka = 1, (b) Ko = 1, Ka = \/2. 




FIG. 9: (color online) Same as Fig. 8, but for an absorbing 
impurity site {Im{Ea) = — 1; other parameters as in Fig. 8) 



eral analysis presented in Sec.II.D, the survival probabil- 
ity P{t) = |ca(i)P decays to zero in the absorbing case 
(Fig. 9), and the existence of a spectral singularity does 
not basically influence the decay dynamics of state \a). 
Conversely, for an amplifying impurity site the occupa- 
tion probability P{t) does not decay to zero. Note that 
for Ha/ 1^0 < \/2 [like in Fig. 8(a)], P{t) converges to a 
steady state value, whereas for Ka/^o — the proba- 
bility P{t) is an oscillating function [see Fig. 8(b)]. The 
different behavior in the two cases is explained by observ- 
ing that, for Kq/ko ^ there is one spectral singularity 
and thus in Eq. (45) there is only one pole that contributes 
to the non-decaying part of Ca{t). On the other hand, for 
the somehow special case Kq/kq = the Haniiltonian H 
has two spectral singularities [see Eq.(57)], and in Eq.(45) 
there are two poles that contribute to the non-decaying 
part of Ca {t) ■ The interference of these two non-decaying 
terms explains the oscillatory behavior of P{t) in Fig. 8(b) 
(see also [35|). 



IV. CONCLUSIONS 



In this work a comprehensive analysis of the spectral 
properties of a non-Hermitian extension of the Friedrichs- 
Fano-Anderson model has been presented. The FFA 
model generally describes the decay of a discrete state \a) 
of energy Ea coupled to a continuum of states. Here we 
have extended the ordinary model by allowing the energy 
Ea to become complex- valued, with either lm{Ea) > 
(the 'amplifying' case) or lm{Ea) < (the 'absorbing' 
case). Contrary to the Hermitian FFA, it has been shown 
by a direct diagonalization procedure and by the analy- 
sis of the resolvent operator that spectral singularities 
in the continuous spectrum may exist for both the am- 
plifying and the absorbing non-Hermitian FFA model. 
The physical implications and relevance of spectral singu- 
larities have been discussed, in particular with reference 
to a tight-binding realization of the non-Hermitian FFA 
model that describes single-particle electronic or pho- 
tonic transport in a semi-infinite lattice attached to an 
impurity site with a complex energy. Different behaviors 
have been found for an amplifying and for an absorb- 
ing impurity site, reflecting the circumstance that the 
divergence of the resolvent G{z) appears when the spec- 
tral singularity is approached either from above or from 
below of the complex energy plane. For an amplifying 
impurity site, the spectral singularity manifests itself as 
a divergence of the reflection probability from the lattice 
boundary, a result which is the 'discrete' analogous of the 
general result recently established by Mostafazadeh for 
wave scattering by complex potentials in the framework 
of the continuous Schrodinger equation [ll|. As com- 
pared to Ref.[ll|, here we have also clarified the phys- 
ical relevance of spectral singularities in the temporal 
domain, showing that in the amplifying non-Hermitian 
FFA Hamiltonian the spectral singularity of the resol- 
vent prevents the decay of state \a) into the continuum, 
in spite of the absence of bound states. For an absorbing 
impurity site, we have shown that a spectral singularity 
corresponds to a zero of the reflection probability from 
the lattice boundary. This result is clearly not observable 
in the problem of wave scattering from complex barriers 
addressed in Refs. [HI, Ull, where the double-degeneracy 
of energy levels plays a major role and a spectral singu- 
larity can never correspond to the total absorption of the 
incident wave. 

Owing to the importance of the FFA model in different 
areas of physics, it is envisaged that the present anal- 
ysis may stimulate further theoretical and experimental 
studies aimed to investigate the unique features of non- 
Hermitian physical systems. In particular, recent results 
obtained in photonic systems based on coupled waveg- 
uides or arrays of coupled waveguides with controlled re- 
gions of optical gain and/or loss [5l|, indicate that en- 
gineered photonic systems might provide an accessible 
laboratory tool to experimentally observe spectral singu- 
larities. 
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Appendix A: Matrix elements of the resolvent 

In this Appendix we derive the expressions of the ma- 
trix elements ga,a{z) = {a\G{z)a), GkMiz) = {k\G{z)a), 
Ga,k{z) = {a\G{z)k) and Gk,k'{z) = {k\G{z)k') of the re- 
solvent [Eqs.(6-9) given in the text]. 
From the identity G(z)(z-iJo-^) = {z~Ho~V)G{z) = 
I it follows that 

{a\G{z){z- Hr,)a) - {a\G{z)Va) = 1 (Al) 

{a\G{z){z - Ho)k) - {a\G{z)Vk) = (A2) 

{k\{z - Ho)G{z)a) - {k\VG{z)a) ^ (A3) 

{k\G{z){z- Ho)k') - {k\G{z)Vk/) = 6(k-k')iA4) 

Taking into account that 



V\a) = / dkv*{k)\k) , V\k) = v{k)\a) 



(A5) 



and that (z — Ho)\a) = (z — Ea)\a), {z — Ho)\k) — {z 
E{k))\k), Eqs.(Al) and (A2) take the form 



{Z - Ea)ga,a (z) - J dkv* {k)ga,k = 1 (A6) 

[z - E{k))ga,k{z) - v{k)gaAz) = (A7) 

which can be solved for and 5a, fc, yielding Eqs.(6) 
and (7) given in the text. To calculate Qk,a{z)^ we 
use Eq.(A3) and note that (fc|(z - HQ)G{z)a) = ((z* - 
Hl)k\G{z)a) - (z - E{k))gk,a{z) and {k\VG{z)a) = 
{Vk\G{z)a) ^ v*ik)ga,aiz)- This yields (z - 
E{k))gk,a{z) — v*{k)ga,a{z) = 0, which can be solved for 
gk,a{z), yielding Eq.(8) given in the text. Finally, the ma- 
trix element gk,k' (z) is obtained from Eq.(A4), which can 
be written in the form {z — E{k'))gk k' {z) — v{k')gk aiz) = 
6{k-k'),i.e. 

. ,._ v{k')gk.a{z) , Sjk^k') 

^'^''^'^ - z-E{k') + J^W)- ^ ^ 

Substitution of Eq.(8) into Eq.(A8) finally yields Eq.(9) 
given in the text. 



Appendix B: Conditions for a real- valued energy 
spectrum of the non-Hermitian Hamiltonian 

In this Appendix we derive the necessary and sufficient 
conditions that ensure a real-valued energy spectrum 
for the non-Hermitian FFA Hamiltonian H introduced 
in Sec. HI. A. As shown in Sec. II. B, this condition is 
equivalent to the vanishing of the point-spectrum of 
H , i.e. to the absence of bound states. The detailed 
calculations can be performed following two different, 
though equivalent, approaches. The first one starts from 
the representation of H in the {|a), |fc)} basis (the Bloch 
basis), whereas the second approach uses a different 
decomposition of H, namely on the {|a), |n)} basis, 



where \n) are the Wannier states introduced in Sec. HI. B 
(the Wannier basis). For the sake of completeness, we 
present the detailed calculations for both approaches. 

1. Bloch-basis representation of H. As shown in 
Sec. II. B, the absence of bound states of H requires that 
Eq.(18) does not admit of any solution in the complex z 
plane. Using the expression (50) of the self-energy S(z), 
Eq.(18) takes the form 



1- I ^ 



Er, 



44 



(Bl) 



We can solve Eq.(Bl) by introducing, in place of z, the 
new complex- valued variable fi defined by 

z = — Ko[exp(/i) -1- exp(— /i)] = ~2ko cosh/j,. (B2) 

Without loss of generality, we may assume Re(^) > 0. 
In fact, the function z(/i) defined by Eq.B(2) is invariant 
for the inversion /i — >■ — /i, so that we may restrict our 
analysis to the case Re(/i) > 0. With such a substitution, 
the square root on the right hand side of Eq.(Bl) can 
be solved analytically, yielding ±2iKQ sinh fi. Some care 
should be taken when choosing the right determination 
(i.e. sign) of the square root [4a]. For Re(/^) > 0, one 
obtains 



2 1- 



2kq 



E 

cosh n H = 1 sinh /i. (B3) 



After setting ^ = exp(/i), from Eq.(B3) one obtains 
Eq.(52) given in the text once cosh/u and sinh/x are 
expressed in terms of the exponentials exp(±/i) = 
Therefore, if the two roots ^1,2 of Eq.(52) satisfy 
the condition 1^1,2] < 1, Eq.(B3) does not have roots 
with Re(/x) > 0, and hence H does not have bound states. 

2. Wannier-basis representation of H. In this ap- 
proach, we use the tight-binding representation of 
the Hamiltonian H using the Wannier function basis 
[Eq.(63)]. Bound states of H correspond in this case to 
surface states localized near the edge of the truncated 
lattice of Fig. 3. They can be directly determined by 
looking for a solution to Eqs. (65-67) of the form 

c„{t) = exp[-^(n-l)-iM], Ca{t) = Aex.p{~iEt) (B4) 

(n > 1), where E is the energy of the surface state. The 
constants /z and A, as well as the dependence of E on 
/i, are readily determined by substituting Eq.(B4) into 
Eqs. (65-67). One obtains 



E = — 2Kocosh^ 
E ~ —Koe:sq){—^) 

EA = -Ka+EaA. 



KaA 



(B5) 
(B6) 
(B7) 



from which the following second-order algebraic equa- 
tions for f = exp(/i) is readily obtained 



(B8) 
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which is Eq.(52) given in the text. Localization of the roots ^i_2 of Eq.(B8) satisfy the condition |^i,2| ^ 1, there 
surface state at the lattice edge requires c„ — as n — ?> are not surface states at the lattice edge, 
oo, i.e. Re(/z) > [see Eq.(B4)]. Therefore, if the two 
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